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ABSTRACT 

An anisotropic viscoplasticity model is developed as an extension of the well 
known Bodner model. The extension is made by replacing the effective stress of the 
isotropic Bodner model by one involving invariants for transverse isotropy. The resulting 
model retains the simplicity of Bodner's in the ease with which the material constants are 
determined experimentally. It allows a representation of strong initial anisotropy yet is 
based on a scalar state variable under the assertion that induced anisotropy is negligible 
relative to the strong initial anisotropy. Temperature dependence is taken as in the original 
Bodner theory. Account is made of fiber volume fraction through nonlinear rules of 
mixture applied to the stress history and anisotropy parameters. Focus is on the theoretical 
development of the model, however, application to a W/Cu composite is in progress and 
will be reported as a sequel to this report. 

INTRODUCTION 

Bodner's model of viscoplasticity [ 1,2 ] is one of several unified models that has 
received attention over the past decade. It was one of two models that was focused on 
under the extensive HOST program sponsored by NASA. In its earliest form, Bodner's 
model applies to initially isotropic materials and incorporates a single scalar state variable. 
In the terminology of continuum mechanics, such a model is labeled isotropically 
hardening. A major strength of the early Bodner model is its simplicity and the relative 
ease with which the material parameters can be determined through experiment. 

The objective here is to develop a simple and tractable transversely isotropic 
model, applicable to metallic composites, by extending the well established Bodner model. 
This is accomplished, in part, by replacing the effective stress ^3 J 2 upon which the 
isotropic Bodner model is based with another effective stress defined in terms of invariants 
reflecting local transverse isotropy. Identification of an appropriate set of invariants for 
transverse isotropy was made by Robinson and Duffy [ 3 ]. The resulting model, capable 
of representing strong initial anisotropy, retains its dependence on a single scalar state 
variable. It is argued that induced anisotropy (e g., through directional or kinematic 
hardening) is negligible relative to the strong initial anisotropy. 



The anisotropic, viscoplasticity model developed here, like the isotropic Bodner 
model, is effective in representing rate-sensitive, non isothermal plastic responses typical 
of the histograms of rocket engines, the space-shuttle main engine (SSME) being a 
prominent example. Temperature dependence is taken as in the original Bodner model. 
Influence of the fiber/matrix volume fraction p is included by considering the internal 
variable (or stress history parameter) Z and the anisotropy parameter C to functionally 
depend on p. 

Here, emphasis is on the theoretical development of the anisotropic model. 
Several plots analogous to those in [ 2 ] are shown illustrating general features of the 
anisotropic theory. An application of the model is in progress to a W/Cu composite for 
which tensile (and creep) data are available over a range of temperatures and strain-rates 
and two volume fractions. The resulting material constants and parameters for the W/Cu 
composite will be reported in a sequel to this report. 

ISOTHERMAL STATEMENT OF THE MODEL 

We first state the multiaxial form of the transversely isotropic model under 
isothermal conditions at a reference temperature T 0 and for a reference fiber volume 
fraction p a . The flow law is 



and the evolution law is 

Z = m(Z s -Z)W W = T^ fi (2) 

in which e tj are the components of inelastic strain rate, cr is the effective stress 

5 = jx/, +(!-<?)/, +|(l-OA) (3) 

and Z is the internal state variable (stress history parameter). 

The invariants /,,/ 2 and / 3 in (3) are defined as 

l, = T-I 2 T=Dfys„ 
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j are the components of deviatoric stress, D y are components of an orientation tensor, cf. 
[ 3 ], e 0 ,n,m,Z s and Z Q (the initial value of Z) are the material constants of 
viscoplasticity; £ and C, are constants relating to the degree of anisotropy. 

The model expressed in (1),(2),(3) and (4) must be supplemented by special 
provisions allowing for stress reversals under cyclic loading conditions. It is expected that 
in metallic composite materials, as with monolithic alloys, a reversal of the sign of stress 
after unloading from an inelastic state is generally accompanied by a marked decrease in 
hardening. In the absence of detailed information on the inelastic cyclic response of 
metallic composites, we shall adopt the same provisions for stress reversals as are outlined 
in [ 1,2 ] for monolithic alloys. 

EQUAL LONGITUDINAL AND TRANSVERSE SHEAR £ = 0 

A simpler form of the model results when the inelastic response in longitudinal and 
transverse shear can be idealized as being equal. This corresponds to taking the material 
parameter £ = 0 . In this case (1) and (2) remain the same, however, (3) and (4) become 



s=^3(y ! -|<7 ! ) 

(5) 

and 


(6) 


respectively. 

We shall use this simplified form in the application to W/Cu. It is interesting to note that 
the effective stress (5) relates directly to the anisotropic yield function used by Robinson 
and Pastor [ 4 ] in the limit analysis of a metallic composite ring. 

FULL ISOTROPY f = 0 AND £ = 0 

Additionally taking £ = 0 the model reduces to the fully isotropic form of Bodner 
[1,2]; (5) and (6) reduce to 

a = y[3^ (7) 
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and 


(8) 


r » = S >J 


As the characterization of the Bodner model is based on uniaxial testing, we 
consider the uniaxial form of the isotropic model represented by (1),(2), (7) and (8). Thus, 



Z = m(Z s -Z)ae (10) 


The constants 

e 0 ,n,Z s ,m and Z a , (11) 

are found by correlating calculated responses based on (9) and (10) with tensile data 
obtained at strain-rates and temperatures of interest. Details of the established 
characterization procedure for the isotropic Bodner model are discussed in [2], 

UNIAXIAL: LONGITUDINAL (0°) STRESS 

We return to the simplified (4 = 0) anisotropic form (1),(2),(5) and (6). Reduction 
to uniaxial stress along the longitudinal (0°) fiber direction gives 



i [ ifrV"' 

— = ex p 0 
£ o 2 \°J 

(12) 


r=m\Z]-Z')crs 

(13) 

with constants 

K> n ,Z*,m* and Z* 

(14) 


We note that (12) and (13) with the constants (14) are identical in form to (9) and 
(10) with the constants (11). Thus, we can use the same characterization procedure on the 
anisotropic material under 0° testing to get the constants (14) as is used on the isotropic 
Bodner model to get the constants (1 1). As indicated earlier, this procedure is well 
established [2], 
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The constants (14) for 0° stress relate to those of the multiaxial model (1),(2),(5) 
and (6) as follows: 





m -m( l~0 



(15) 


. , _ <J , A o 

Fig. 1 shows the dependence of the uniaxial flow stress parameter = — — — under 0 


stress on the strain rate parameter — = ■ for different values of the constant n . 

With £ = 0 (isotropy) Fig. 1 reduces to that given in [ 2 ]. 


Thus, through (15), the multiaxial model is totally specified once the constants 
(14) have been determined by 0° testing and the value of the anisotropy parameter C, is 
known. Determination of C, will be discussed in a later section. 


UNIAXIAL: TRANSVERSE (90°) STRESS 

Next, we specialize the (4 = 0) model (1),(2),(5) and (6) for transverse (90°) 
stress.There results 


e 

— =exp 



Z~=/n~(Zr-Z~)(Tf 

(H) 

with constants 

e~ ,n,Z~ ,m~ and Z~ 

(18) 
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(16) 


Again, the form of (16),(17) with (18) is identical to that of the uniaxial Bodner model 
(9),(10) with (1 1), so the same characterization procedure also can be applied to 90° 
testing. The resulting constants (18) for 90° relate to the constants in the multiaxial model 
(1),(2),(5) and (6) according to: 
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( 19 ) 


A plot identical to Fig. 1 but for 90° stress is obtained by replacing by and by 

z z s 0 

— — . Again, the multiaxial model is totally specified once the constants (18) are 
determined from 90° testing and the value of C, is known. 

DETERMINATION OF £ 

Comparing (15) and (19) that relate the constants for 0° and 90° testing to those 
of the multiaxial model, it is apparent that £ is specified once the sets of constants (14) 
and (18) are known. However, both sets need not be complete. If, for example, we know 
all the constants (14) from 0° testing and have only, say Z~ from 90° testing, C, can be 
determined from 



( 20 ) 


We select Z~ (the initial value of the internal variable Z~) in this context because it is one 
of the most easily obtained constants in the characterization procedure. The same could be 
done using the constant s ", in which case C, would be obtained from 

r? 

— 1 ( 21 ) 

l-c 



We note from (20) and (21) that 
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Evidently, C, is over specified in that four equations like (20) and (21) can be 
written from (15) and (19). It is possible that one value of C, cannot be found satisfying all 
four equations to a desired degree of accuracy. Our recourse then is to reexamine some of 
the idealizations made in the simplified (£ = 0) model (1),(2), (5) and (6). An obvious 
consideration is to remove the condition £—0 leaving the degree of anisotropy to depend 
on both £ and C' the model is then (1),(2),(3) and (4), as originally stated. The simple 
characterization procedure just outlined, based directly on that established for the Bodner 
model, is not generally applicable to the more comprehensive model. An analogous but 
more involved procedure has been identified for this model; unavoidably it requires 
transverse and longitudinal shear testing to determine £ . As this kind of test data on 
metallic composites, in particular W/Cu, are not readily available at present, we shall 
continue to assume £ = 0 and to develop and apply the simplified model (1),(2),(5) and 
(6). Development and characterization of the more general anisotropic model including 
both £ and C, will be addressed as a continuation of this research. 

In the context of the simplified (£ = 0) model. Figs. 2-4 show the influence of the 
anisotropy parameter C on the transverse (90°) stress response. Although these figures 
relate to transverse behavior, they are shown in terms of the uniaxial flow stress parameter 


— and the strain rate parameter -4- for comparison with the 0° response. Figs. 2,3 and 
Z* K 

4 correspond to £= 0.3, 0.6 and 0.9, respectively. We see that the 90° flow stress 
parameter decreases as C, increases for fixed values of strain rate and n . In particular, for 


C = 0.9 and n = 2 (Fig. 4), * 0.2 when 4- = 10 

Z e„ 


,-2 


Comparing this with longitudinal 


(0°) stressing for n = 2 (Fig. 1), * 0.6 when = 10‘ 2 . Thus, for the degree of 

Z £ 0 

anisotropy £ = 0.9 and with n = 2, the ratio of the flow stress for 0° to that for 90° is 

S _2 

a 3 for an inelastic strain rate of — = 10 . 
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TEMPERATURE DEPENDENCE 


As with the isotropic Bodner model, the primary influence of temperature on the 

C 

flow stress is taken through the parameter n(T ) . If we take n{T) = — where k is 

fv A 

Boltzmann's constant, a thermal activation form results with the anisotropic activation 
energy function 




(23) 


With full isotropy, 4 = 4 = 0, (3) or (5) reduce to (7), and (23) becomes the activation 
energy function specified in [ 2 ]. 

Fig. 5 shows the dependence of the flow stress parameter for 0° uniaxial stress on 


— = for different values of the strain rate parameter. Again, for 4” = 0 (isotropy) this 

n C! k 

plot reduces to its counterpart given in [ 2 ]. . 

It may prove necessary to include temperature dependence in other material 
parameters, eg., Z a ( the initial value of the internal variable). However, in the application 
to W/Cu we shall attempt to limit temperature dependence to the parameter n . 


FIBER VOLUME FRACTION p 


The characterization procedure described above is defined as being applicable to a 
composite material having the reference fiber volume fraction p a . A composite having the 
same constituents but a different fiber volume fraction p can be viewed as a distinct 
material and the same characterization procedure applied to it, yielding another set of 
constants and parameters. In principle, this can be repeated for a range of fiber volume 
fractions of interest with the resulting constants fit in some optimal sense as functions of 
p. This may be necessary in practice because the precise fiber volume fraction is often not 
known a priori for a structure fabricated from a composite material. Indeed, p may vary 
from point to point in a given body or structure. A constitutive model to be used in the 
analysis of such structures must allow for variations in p. 

To avoid a large amount of tedious and costly characterization testing, we can 
appeal either to micromechanical studies to provide guidance of how the overall 
composite response changes with fiber volume fraction p, or we can attempt to limit the 
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amount of phenomenological testing by judiciously identifying the key parameters in which 
the p dependence should reside. Here, we pursue the latter course within the context of 

the simplified £ = 0 model (1),(2),(5) and (6). 

Observations based on data for two fiber volume fractions of a W/Cu composite 
suggest that the primary dependence on p should be taken in the state variable (or stress 
history parameter) Z*and in the anisotropy parameter These dependencies are taken as: 

Z\p) = Z'{p 0 )'¥(p) (24) 

ap)=ap 0 )®(p) ( 25 > 


where p Q is the reference fiber volume fraction and and O can be considered 
nonlinear rules of mixture. Of course, Z\p 0 ) in (24) is identical to Z* in (12)-(14). 

Consistent with the basic features of the model, exponential forms are chosen for 
*¥ and O as follows: 


V(P) = 


1 - 


Z\0) ' 
Z\p 0 ) J 


1 - exp( — 


ap 


Po~P ) 


Z\ 0) 
Z'{p 0 ) 


(26) 


O(p) = 1 - exp( — (27) 

Po~P 

in which a and /? are constants. 

From (24)-(27), we see that as p p 0 , Z’ Z\p a ) and £(p 0 ), relating 
as earlier, to the reference volume fraction. As p — > 0 (no fibers), Z — > Z (0) 
corresponding to the matrix material, and £ -» 0 corresponding to isotropy. 


Figs 6 and 7 are plots of '¥(p)vs — and <h(p)vs— , respectively, where for 

Po Pc 

illustration we have taken Z* (0) / Z*(p Q ) = l/ 3, a = l/ 2 and ^ = 3. Fig. 8 shows the 

dependence of the flow stress parameter — ^ on the strain rate parameter -77- for n — 1 and 

Z & o 

2 1 

for three values of the fiber volume fraction p — p 0 , ~^Po “A,- Values of the other 

relevant parameters in Fig. 8 are taken as in Figs. 6 and 7. The three upper curves in Fig. 8 
are for longitudinal (0°) stress; the three lower curves are for transverse (90°) stress. As 
expected intuitively, both the 0° and 90° flow stress parameters decrease with decreasing 
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fiber volume fraction. Also as expected, the change in the flow stress with p is 
considerably greater for 0° than for 90° (« 4 times greater in this example ). 

The additional constants introduced in (24)-(27), viz., Z*(0) (or equivalently, 

Zj(0) and Z*(0)- the initial and saturation values of the state variable for the matrix ), a 
and /?, must be determined using experimental results from tests on, at least, one fiber 
volume fraction other than the reference p a . 

SUMMARY AND CONCLUSIONS 

A simple and tractable, transversely isotropic viscoplasticity model is developed by 
extending the well known isotropic Bodner model. The model retains the simplicity of the 
Bodner theory regarding the relative ease with which the material constants and 
parameters can be determined. It is capable of representing strong initial anisotropy yet is 
based on a single scalar state variable under the assertion that induced anisotropy through 
directional or kinematic hardening can be ignored relative to the strong initial anisotropy. 
Like the original Bodner model, the anisotropic model is particularly effective in 
representing rate-sensitive, non isothermal viscoplastic responses typical of the histograms 
of rocket engines, e.g., the SSME. 

A procedure for determining the constants of the model is discussed and is shown 
to be essentially identical to the well established characterization procedure for the Bodner 
model. Temperature dependence is taken as in the original Bodner theory. Account is 
made of varying fiber volume fraction by including its dependence specifically in the state 
variable Z ( stress history parameter) and the anisotropy parameter C , . 

Application of the model to a W/Cu composite is in progress. The results 
specifying the material constants and parameters will be reported as a sequel to this report. 
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Fig. 1. Dependence of the uniaxial flow stress parameter on the strain rate parameter 
under 0° stress for different values of n . 
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Fig. 2. Dependence of the uniaxial flow stress parameter on the strain rate parameter 
under 90° stress for different values of n and C = 0 3. 



Fig. 3. Dependence of the uniaxial flow stress parameter on the strain rate parameter 
under 90° stress for different values of n and £ = 0.6. 
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Fig. 4. Dependence of the uniaxial flow stress parameter on the strain rate parameter 
under 90° stress for different values of n and £ = 0 . 9 . 



n Clk 


Fig. 5. Dependence of the uniaxial flow s' ress parameter on — = under 0° stress for 

n Clk 


different values of the strain rate parameter. 


13 








Fig. 8. Dependence of the uniaxial flow stress parameter on the strain rate parameter 

under 0° and 90° stress with n-2 and p- p o , —p o and - p o . 

3 3 
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